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ABSTRACT 


VOID  FRACTION  IN  PARALLEL  CHANNEL  BOILING  WATER  REACTORS  WITH 
APPLICATION  TO  'TWO-PHASE  TURBULENT  MIXING 

RONALD  MALCOLM  REESE 


Submitted  to  the  Department  of  Naval  Architecture  and 
Marine  Engineering  on  June  4,  1970,  in  partial  fulfillment 
of  the  requirements  for  the  degrees  of  Naval  Engineer  and 
Master  of  Science  in  Mechanical  Engineering. 


An  experiment  using  air-water  mixtures  to  model  steam- 
water  flow  in  a  reactor  was  devised  to  determine  the  turbulent 
mixing  mass  flow  rates  between  two  subchannels  over  a  wide 
range  of  flow  conditions  and  therefrom  predict  turbulent 
mixing  occurring  in  the  reactor.   An  experimental  apparatus 
was  constructed  for  this  experiment. 

Void  fraction  was  chosen  as  a  plausible  correlating 
parameter  for  turbulent  mixing  and  was  measured  together  with 
axial  pressure  gradient  for  two  subchannel  geometries  of 
interest . 

The  subchannel  geometries  investigated  were  those  found 
in  a  typical  BWR  fuel  bundle.   The  flow  regime  in  the  sub- 
channels was  detected  by  means  of  an  electric  conductivity 
probe  and  visual  observations.   Data  obtained  are  compared 
with  various  correlations. 

Predictions  of  void  fraction  and  pressure  gradient  with 
the  Lockhart-Martinelli  correlation  agreed  closely  with  the 
data.   Comparisons  were  also  made  with  Griffith's  slug  flow 
model,  and  with  the  Cravarolo-Hassid  and  Marchaterre-Hoglund 
void  fraction  correlations. 
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A  =  subchannel  cross  sectional  area  (FT  ) 

c  =  rod  gap  width  (FT) 

two-phase  velocit 
parameter  (dimensionless) 


C   =  two-phase  velocity  and  concentration  distribution 


-=£■)    .  =  frictional  pressure  gradient  if  the  liquid  phase 
flowing  alone  were  allowed  to  fill  the  entire 
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-~)  __  m_,  =  frictional  two-phase  pressure  gradient  (LB  ../FT  -FT) 
dl  FR  TP  c  sr  a  f 

■v?-)  „_  _„  m_  =  two-phase  elevation  pressure  gradient  (LB  ../FT  -FT) 
dl  ELEV  TP         ^  tr  -a  f 

D   =  hydraulic  diameter  (4A/P)  (FT) 

f  =  single  phase  liquid  friction  factor  (dimensionless) 

f  ,  f   =  diversion  cross  flow/turbulent  mixing  momentum 

exchange  multiplier 

F  =  friction  force  per  unit  length  at  the  subchannel  wall 
(LBf/FT) 
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g  =  acceleration  due  to  gravity  (FT/HR  ) 
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g   =  conversion  factor  (LB  FT/LB..HR  ) 
^c  m   /   f 
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G  =  mass  velocity  or  mass  flux  (LB  /HR  FT  ) 

G.,  G.  =  mass  velocity  in  subchannel  i ,  j  in  the  axial 

direction  (LB  /HR  FT2) 
m 

G. .',  G..1  =  turbulent  mixing  mass  velocity  from  subchannel 

ID     Di  y  *       2 

i  to  j  and  vice  versa  (LB  /HR  FT  ) 
_i  m 
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G  .  .    =  net  turbulent  mixing  mass  velocity  (LB  /HR  FT  ) 
mixing  ?  J     m 

h  =  enthalpy  (BTU/LB  ) 

3       2 
j  =  local  volumetric  flux  density  (FT  /HR  FT  )  =  V  a   +  Vf (1-a) 

3       2 

3      =    local  gas  volumetric  flux  density  (FT  /HR  FT  ) 
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£  =  mixing  length  (FT) 

L  =  communication  (window)  length  (FT) 

m. ,  m.  =  inlet  mass  flow  rate  in  subchannel  i ,  j  or  exit  mass 
i    3 

flow  rate  under  isolated  conditions  (LB  /HR) 

m.*,  m.*  =  exit  mass  flow  rate  from  subchannel  i,  j  when 
subchannels  are  communicating 
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p  =  pressure  (LBf/FT') 
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Q  =  volume  flow  rate  (FT  /HR) 

q  =  heat  flow  per  unit  length  of  subchannel  (BTU/HR  FT) 

S  =  slip  ratio  (V  / Vf ,  local;  V  /Vf ,' average) (dimensionless) 

u  =  axial  velocity  (FT/HR) 

u.  .   ;  u..    =  diversion  cross  flow  velocity  from  subchannel 

i  to  j  and  vice  versa  for  gas  and  liquid 
Uij,f;  Uji,f    respectively  (FT/HR)' 


u.  .     =  representative  measure  of  the  fluctuating  turbulent 
mixing  gas  velocity  from  subchannel  i  to  j  (e.g., 
RMS  value  or  average  of  absolute  value) (FT/HR) 


u.  .  r-\    =  representative  measure  of  the  fluctuating  turbulent 
1  i],f  '      * 

mixing  liquid  velocity  from  subchannel  i  to  j  (e.g., 
RMS  value  or  average  of  absolute  value) (FT/HR) 
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v'  =  mean  f lowing  two-phase  specific  volume  (FT  /LB  ) , 
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V  =  local  liquid  velocity  (FT/HR) 

V  =  local  gas  velocity  (FT/HR) 

-  Q     Q 

V  =  weighted  mean  gas  velocity  (FT/HR)  =  r-2-  =  — ~ 

-  g2f"   °f 

V_  =  weighted  mean  liquid  velocity  (FT/HR)  =  —  =  ■  ,.,   > 

V  .  =  local  gas  drift  velocity  (V  -j) (FT/HR) 

V  '.    =   weighted  mean  gas  drift  velocity  (FT/HR) 


w. .,  w..  =  diversion  cross  flow  mass  flow  rate  per  unit 
ID    Di 

length  from  subchannel  i  to  j  and  vice  versa 

(LB  /HR  FT) 
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.  .  ' ,  w . .  '  =  turbulent  mixing  mass  flow  rate  per  unit  length 

from  subchannel  i  to  j  and  vice  versa  (LB  /HR  FT) 

J  m 

w1    =  net  turbulent  mixing  mass  flow  rate  per  unit  length 

(LB  /HR  FT) 
nr 

X  =  flowing  quality,  (m  /m  +mf ) (dimensionless) 

a  =  void  or  gas  volumetric  fraction,  (A  /A) (dimensionless) 

3  =  average  mass  velocity  multiplier  for  turbulent  mixing  in 
the  COBRA  Code  (Mixing  Stanton  Number) (dimensionless) 

3 

p  =  average  density  (LB  /FT  )  (pf(l-a)  +  p  a) 

2    2     . 
(f)  „  ,  6=  liquid  and  gas  two-phase  frictional  pressure 
x,    g 

multipliers  (dimensionless) 
<  >  =  denote  cross  sectional  averages,  <f>  =  —  /  fdA 

Subscripts 

eq  =  refers  to  equilibrium  conditions 

f  =  refers  to  liquid  phase 

g  =  refers  to  gas  phase 

i  =  refers  to  ith  subchannel 

j  =  refers  to  jth  subchannel 
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CHAPTER  I 


Introduction 


The  core  of  modern  boiling  water  reactors  is  composed  of 
a  large  number  of  subassemblies  or  fuel  bundles.   Each  sub- 
assembly contains  a  matrix  of  vertical,  cylindrical,  parallel 
fuel  rods.   The  fuel  bundles  are  generally  enclosed  in  a 
surrounding  unheated  wall.   Figure -1  shows  such  a  typical 
16-rod  bundle  arranged  on  a  square  array.   The  main  coolant 
flow  is  vertical,  parallel  to  the  axis  of  the  subassembly. 
However,  secondary  radial  flows  between  adjacent  rods  and 
between  the  outer  rods  and  the  wall  also  take  place. 

The  assumption  is  generally  made  that  the  phenomenon 
known  as  critical  heat  flux  (CHF) ,  burnout,  or  departure  from 
nucleate  boiling  (DNB)  depends  only  upon  the  local  conditions. 
The  prediction,  then,  of  the  location  and  the  magnitude  of  the 
CHF  requires  a  complete  and  detailed  knowledge  of  the  main  and 
secondary  flow  patterns. 

Accurate  estimation  of  the  magnitude  of  the  secondary 
flows,  known  as  mixing,  has  been  the  major  problem  in  the 
prediction  of  the  critical  heat  flux  condition.   Considerable 
progress  has  been  made  in  the  area  of  mixing  flows  in  single 
phase  flow.   The  mechanisms  of  two-phase  mixing,  however,  are 
still  not  well  understood. 

It  has  been  convenient  to  subdivide  the  main  flow  through 
the  rod  bundle  or  subassembly  into  flows  between  four  adjacent 
rods  or  between  a  rod  or  rods  and  the  subassembly  wall.   Such 
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artificially  defined  flow  ducts  are  generally  referred  to  as 
rod  channels  or  subchannels. 

The  standard  approach  to  a  thermal-hydraulic  analysis 
of  a  rod  bundle  utilizes  these  subchannels.   A  typical 
subchannel  breakdown  of  a  rod  bundle  is  illustrated  in  Figure 
1.   Mean  properties  are  assigned  to  each  subchannel.   Thus  a 
mean  axial  velocity,  enthalpy,  void  (or  volumetric)  fraction 
is  assumed  to  exist  over  the  entire  cross  section  of  a  sub- 
channel.  This  lumped  parameter  approach,  although  far  from 
being  strictly  correct,  serves  as  a  first  approximation  to 
the  actual  velocity,  enthalpy  and  concentration  profiles  ex- 
isting in  a  given  subchannel. 

Mixing  between  parallel,  ventilated  (communicating)  rod 
channels  has  been  conveniently  defined  by  Todreas  and  Rogers 
(2)  as  follows: 

Natural  mixing 

a)  Turbulent  interchange 

b)  Diversion  cross  flow 
Forced  mixing 

a)  Flow  scattering 

b)  Flow  sweeping 

Turbulent  interchange  is  caused  by  the  fluctuating 
velocities  and  pressures  present  in  turbulent  flow  and  has 
been  broadly  defined  as  a  no  net  mass  transfer  phenomenon,  in 
accordance  with  Prandtl ' s  mixing  length  theory.   Such  a 
definition  has  application  to  single  phase  flow  but  requires 
revision  for  two- phase  flow. 

11  • 
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Diversion  cross  flow  is  the  net  secondary  flow  resulting 
from  a  static  pressure  difference  between  parallel  channels  at 
a  given  level  within  the  rod  bundle.   In  general,  it  will  be 
represented  by  a  net  mass  transfer  for  single  and  two-phase 
flow  conditions. 

Forced  mixing,  by  flow  scattering  and  flow  sweeping,  is 
artificially  induced.   Rod  spacers,  for  example,  may  produce 
a  random  diversion  of  flow  (flow  scattering) .   Devices  such  as 
helical  rod  wraps  produce  a  preferential  diversion  of  flow 
(flow  sweeping) . 

In  the  realm  of  natural  mixing  it  has  been  assumed  that 
the  mechanisms  of  turbulent  interchange  and  diversion  cross 
flow  are  independent  or  at  least  may  be  analyzed  as  such  as  a 
first  approximation.   As  noted  by  Todreas  and  Rogers  (2), 
where  large  diversion  cross  flows  exist,  this  probably  is  an 
erroneous  assumption,  but  nevertheless,  it  has  been  a 
reasonable  method  of  attacking  the  mixing  problem. 

The  magnitude  of  the  mixing  flows  has  been  investigated 
by  Rowe  and  Angle  (1) ,  (8)  and  by  Lahey ,  Shiralkar  and 
Radcliffe  (3)  for  the  diabatic  case.   Rowe  and  Angle 
investigated  turbulent  interchange  by  using  radioactive  tracers 
This  method  could  not,  however ,  give  the  net  mixing  vapor  and 
liquid  mass  flow  rates  between  two  channels.   It  has  been 
shown  (e.g.  (3))  that  the  vapor  has  an  affinity  to  migrate  to 
the  less  restricted  regions  of  the  central  channels  away  from 
the  unheated  bundle  walls. 
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It  v/ould  appear,  then,  that  two  phase  natural  mixing  in 
the  absence  of  a  radial  pressure  gradient  may  be  caused  by  at 
least  two  phenomena:   1)  turbulent  interchange,  as  in  single- 
phase  flow,  and  2)  a  bubble  drift  to. less  restricted  flow 
areas . 

In  spite  of  this  information,  in  this  thesis  natural 
mixing  flow  in  the  absence  of  a  radial  pressure  gradient  will 
continue  to  be  referred  to  as  turbulent  mixing.   This  is  in 
concurrence  with  most  literature. 

The  mechanism  of  diversion  cross  flow  seems  to  be  fairly 
well  understood,  and  the  magnitudes  of  these  flows  can  be 
estimated  by  simply  modeling  the  rod  gap  as  an  orifice  subjected 
to  a  pressure  difference. 

The  mechanisms  of  turbulent  mixing  flows,  on  the  other 

hand,  are  not  at  all  well  understood,  and  past  estimates  have 

often  been  based  on  crude  and  somewhat  arbitrary  models.   A 

typical  method  has  estimated  the  turbulent  mixing  mass  velocity, 

G  .  .   ,  as  equal  to  some  constant,  3,  multiplied  by  the 
mixing       -1  '   '       l  J 

average  mass  velocity  of  adjacent  subchannels. 

w! .     (G.  +  G.) 

G  .  .    =  -il  =  g— L-— JL- 

mixing     c  2 

where  w! .  is  the  mixing  mass  flow  rate  per  unit  length  and  c  is 
the  rod  gap  spacing. 

With  this  background,  it  can  be  seen  that  a  more  accurate 
determination  of  the  magnitude  and  mechanisms  of  turbulent 
mixing  flows  is  required. 
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It  was  decided,  then,  to  construct  an  apparatus  which 
could  determine  the  magnitude  of  the  turbulent  mixing  flows 
over  a  wide  range  of  flow  conditions.   This  apparatus, 
hopefully,  could  also  give  some  insight  into  the  mechanisms 
of  two -phase  turbulent  mixing.   Such  an  apparatus  would  have 
to  exclude  the  possibility  of  diversion  cross  flow  mixing 
and  thus  measure  only  the  turbulent  mixing  component.   This 
could  be  accomplished  if  pressure  level  and  pressure  drop 
could  be  equated  at  corresponding  points  in  two  adjacent 
subchannels,  thus  removing  the  driving  force  for  diversion 
cross  flow.   As  a  first  approximation,  such  an  experiment 
could  be  conducted  under  adiabatic  conditions  with  air  and 
water,  modeling  the  steam-water  flow  of  the  reactor. 
Particular  subchannels  of  interest  are  shown  in  Figure  1. 
Subchannel  geometry  (b)  was  selected  initially  and  is  that 
geometry  investigated  in  this  thesis. 

It  was  furthermore  postulated  that  void  fraction  might 
serve  well  as  a  correlating  parameter  for  turbulent  mixing. 
An  indication  of  the  influence  of  void  fraction  on  mixing 
is  the  experimental  evidence  obtained  by  Rowe  and  Angle  (8) 
and  Lahey  et  al.  (3) ,  who  showed  mixing  to  be  a  function  of 
quality,  peaking  somewhere  in  the  slug  flow  regime. 

The  purpose  of  this  research  work  was  the  experimental 
determination  of  void  fraction  over  a  wide  range  of  air  and 
water  mass  flow  rates  for  a  typical  rod  bundle  subchannel 
geometry.   Simultaneously,  pressure  level  and  pressure  drop 
were  recorded  and  flow  regimes  were  identified.   Thus  a 
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complete  picture  of  the  flow  conditions  was  determined  for 
a  wide  range  of  operating  conditions.   The  data  obtained 
have  been  compared  with  various  correlations ,  and 
similarities  and  differences  are  discussed. 


16 


CHAPTER  II 


Theom 


It  is  generally  accepted  that  steam-water  mixtures  may, 
under  certain  conditions,  be  accurately  modeled  by  air-water 
mixtures.   In  the  present  case  this  greatly  simplifies  the 
equipment  necessary  for  the  study  of  two-phase  turbulent 
interchange. 

As  previously  stated,  turbulent  interchange  is  defined  as 
the  mass  transfer  occurring  in  the  absence  of  a  static  average 
pressure  difference  between  adjacent  subchannels  at  a  given 
axial  position.   Hence,  any  experiment  designed  to  determine 
this  mixing  flow  would  necessarily  require  equating  pressure 
level  and  pressure  drop  at  corresponding  points  in  the  sub- 
channels.  By  establishing  flows  meeting  these  conditions  with 
unequal  void  fractions  in  adjacent  subchannels,  one  could 
observe  the  net  mixing  mass  flow  rates  of  air  and  water. 

The  mixing  flow  rates  could  be  obtained  if  some  means 
were  provided  for  isolating  the  two  subchannels  upstream  and 
downstream  of  a  communication  length,  by  simply  measuring  the 
flow  rates  of  air  and  water  in  the  two  subchannels .   This 
method  is  probably  most  accurate  when  the  flows  in  the 
adjacent  subchannels  have  distinctly  different  qualities. 

The  net  mixing  mass  flow  rates  could  also  be  determined 
from  the  difference  in  the  inlet  and  exit  air  and  water  mass 
flow  rates. of  each  channel,  with  the  channels  continuously  in 
communication.   However,  the  maintenance  of  unequal  void 
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fractions  and  yet  equal  pressure  and  pressure  gradients  under 
these  conditions  is  difficult  if  not  impossible.   Furthermore, 
in  the  first  method,  conditions  of  fully  developed  flow  can 
be  attained  prior  to  allowing  mixing  to  occur.   This  was 
thought  to  be  a  more  realistic  approach. 

The  operating  procedure  then  is  to  completely  isolate  the 
subchannels,  establish  flows  in  each  subchannel  satisfying  the 
conditions  of  equal  pressure  and  pressure  gradient,  and  measure 
the  flow  rates.   Subsequently,  a  portion  of  the  common  boundary, 
the  rod  gap,  is  opened  (that  is,  the  subchannels  are  allowed 
to  communicate  over  a  variable  length,  the  "window").   With 
the  "window"  open, the  exit  flow  rates  are  measured  again  and 
the  net  mixing  flows  of  air  and  water  are  calculated. 

This   net  liquid  mixing  flow  is  then  given  by: 

w'  ,  _  *  L  =  (wl .  _-w! .  £)*L  =  m  *  -m.   =-(m  *  -m.  _) 
net,f         31, f   13, f        i,f   i,f     ;j,f   j,f 


where 


1  .  j.   =  net  liquid  turbulent  mixinq  mass  flow  rate  per 
net,f  a  3  c 

unit  length  of  channel  (LB  /HR  FT) 


w 
ne 


w! .  f   =  turbulent  mixing  liquid  mass  flow  rate  per  unit 

length  from  subchannel  i  to  subchannel  j 

(LB  /HR  FT) 
m 


ji,f 


turbulent  mixing  liquid  mass  flow  rate  per  unit 

length  from  subchannel  j  to  subchannel  i 

(LB  /HR  FT) 
m 


L       =  length  of  communication  (FT) 

m.  j.;'  m.  ,  =  exit  mass  flow  rates  of  liquid  from  the 

ith  and  jth  subchannels  under  isolated 

conditions  (LB  /HR) 
m 
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m.  *fi    m-*f   =■   same  as  m.  ,.;  m .  _,  except  that  subchannels 
are  allowed  to  communicate 

A  similar  expression " can  be  written  for  the  net  gas 

turbulent  mixing  mass  flow  rate. 


*  _ 


m 


w^-i-  rr    '    L  =  (w'n         -  w!  .   )  -L  =  m.  *  -m.   =-  (m. 

net,g  HI, g      13 ,g  i,g   i,g      3,9   Drg 

The  net  total  mixing  mass  flow  rate  •  is  w1    =  w*    _  +  w1 

net     net,f     net,g 

The  window  length,  L,  will  necessarily  be  chosen  through 
experimentation.   Hopefully,  by  measuring  the  mixing  mass  flow 
rates  with  decreasing  lengths  of  communication,  differential 
values  can  be  obtained  at  the  limit  by  extrapolation  to  a  zero 
window  opening. 

Subchannel  geometry  (b)  Figure'  1,  selected  for  initial 
investigation,  is  representative  of  a  typical  wall  subchannel 
on  the  left  side  of  the  rod  gap.   The  right  side  of  the  rod 
gap  models  approximately  a  central  rod  bundle  subchannel,  such 
as  that  existing  between  four  rods. 

A  recognized  artificiality  in  this  experiment  is  the 
effect  of  the  imposition  of  a  wall  at  the  rod  gap  to  isolate 
the  channels.   Since  the  width  of  the  rod  gap  (.168")  is  small 
in  comparison  to  the  total  perimeter  (2.071"  and  2.301"  for 
the  left  and  right  channels,  respectively),  it  may  be  argued 
that  the  effect  on  the  overall  flow  conditions  is  small. 
However,  the  abrupt  termination  of  this  wall  at  the  beginning 
of  the  communication  window  in  the  mixing  apparatus  will  have 
some  effect  upon  mixing .   The  eddies  produced  by  the  vorticity 
in  the  boundary  layer  undoubtedly  will  cause  some  additional 


20 


mixing.   It  is  argued,  however,  that  with  sufficiently  long 
window  opening,  the  true  turbulent  mixing  flows  to  be  measured 
will  dominate  the  artificially  induced  effects  due  to  boundary 
layer  vorticity. 

A  further  effect  of  the  presence  of  an  isolating  wall 
between  the  subchannels  is  the  deviation  of  the  void  fraction 
distribution  from  that  which  would  exist  if  no  wall  were 
present  at  the  rod  gap.   It  is  expected,  also  because  of  the 
narrow  rod  gap,  that  this  effect  will  be  small.   Nevertheless, 
if  the  two-phase  mixing  length  is  found  to  be  on  the  order  of 
magnitude  of  the  rod  gap  width,  (i.e.,  if  the  average  sub- 
channel properties  are  not  representative  of  the  local  conditions 
where  mixing  occurs)  concentration  profiles  in  the  communicating 
channels  will  probably  have  to  be  investigated. 

The  net  mixing  mass  flow  rate  per  unit  length  for  gas  and 
liquid  combined  can  be  postulated  to  be  a  function  of  the 
average  void  fraction,  a,'  as  follows: 

w'  ,  =  w!.-w!.  =c[(a.-a.)pjr|u..  ,.  1  +  (a  .-a  .  )  p  lu .  .   I] 
net     31   i]        1   j  Kf  13,  f1    j      1  yg'  13, g1 

In  this  formulation  c  is  the  rod  gap  spacing  and  |u. .  f |  and 

13,1 

lu. .     are  representative  measures  of  the  turbulent 
fluctuating  velocities  of  the  liquid  and  gas  respectively, 
e.g.,  RMS  values  or  averages  of  the  absolute  values. 

The  conservation  equations  and  the  turbulent  mixing 
equation  given  above  are  developed  in  detail  in  Appendix  A. 
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Since  the  net  liquid  and  gas  mixing  flows  are  expected  to 
be  some  function  of  the  void  fractions  in  the  communicating 
channels,  a  detailed  knowledge  of  void  fraction  is  required 
over  a  wide  range  of  water  and  air  mass  flow  rates.   In  order 
to  match  pressure  level  and  pressure  drop  at  unequal  void 
fractions,  a  knowledge  of  pressure  drop  as  a  function  of  void 
fraction  would  also  be  desirable.   Lockhart  and  Martinelli  (5) 
have  shown  the  void  fraction  and  pressure  gradient  to  be 
functions  of  gas  and  liquid  density,  and  hence,  functions  of 
pressure  level.   For  this  reason,  the  pressure  level  at  the 
location  of  the  communicating  window  is  also  of  interest. 

Void  fraction  data  for  air-water  and  other  two-phase 
mixtures  in  horizontal  round  tubes  has  been  correlated  by 
Lockhart  and  Martinelli  (5)  using  a  parameter,  X,  which  is  a 
function  of  flowing  quality,  density  and  viscosity  of  the  two 
phases.   This  same  parameter,  X,  is  also  used  to  correlate  the 
two-phase  frictional  pressure  drop,  dp/dl  ).,-,„     to  the  single- 

phase  gas  or  liquid  pressure  drop  through  the  pressure  drop 

2 

multiplier,  <j>  ,  where 

dp.       _  ,2  dp.      __  .2  dp. 
di'FR  TP    9£  dl'FR  £    9g  dlJFR  g 

The  subscripts  T.P.,  I,    and  g  refer  to  two-phase,  liquid 
only,  or  gas  only  (respectively)  filling  the  entire  tube. 
Thus,  for  given  liquid  and  gas  mass  flow  rates ,  it  is  possible 
to  predict  the  void  fraction,  a,  and  the  two-phase  frictional 
pressure  drop. 
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The  elevation  pressure  drop  is  given  by 


d£ 


dl}ELEV  TP  =   p9/gc  =  [pf  U~a)+  Pga]  gAjc 

Hence,  using  the  Lockhart-Martinelli  correlation,  the  total 
pressure  drop  may  be  predicted. 


dl' TOTAL  TP    dl'FR  TP    dlyELEV  TP 


Void  fraction  correlations  of  other  authors  (12),  (13)  are  also 
used  in  this  work. 

Another  parameter  of  interest  in  two-phase  flow  is  the 
ratio  of  average  gas  velocity  to  average  liquid  velocity, 
known  as  the  velocity  ratio  or  slip'  ratio. 
V 

V£ 

As  noted  by  Zuber  and  Findlay  (9) ,  in  a  flow  channel,  local 

slip  ratios  exist  as  well  as  a  local  concentration  (void 

fraction)  profile.   Knowledge  of  local  void  fraction  may  be 

helpful  in  the  determination  of  a  mixing  length  for  turbulent 

mixing.   Zuber  and  Findlay  have  defined  a  parameter,  C  , 

which  is  a  measure  of  the  peaking'  of  the  concentration  and 

velocity  profiles.   Average  void  fraction  data  for  a  wide 

range  of  flow  rates  provides  sufficient  information  for 

estimation  of  C  . 

o 
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For  round  tubes,  Zuber  and  Findlay  assumed  velocity  and 
concentration  distributions  as  follows: 


4-  =  1  -  (^-)m 


!-(£>" 


a  -a  R 

c   w 


where  j  is  the  local  volumetric  flux  density  of  the  mixture 

3       2 

(FT  /HR  FT  )  and  the  subscripts  c  and  w  refer  to  conditions  at 

the  centerline  and  the  wall  respectively. 

The  distribution  parameter  C   is  defined  as 

^  o 

1/A  /_ajdA  .  _       a 

c     = A__ =  <a3>    _  i  ,   2   n   w  > 

o    [l/A  /  adA] • [1/A  /AJdA]    <a>  <j>      m+n+2  v    <a>; 

where  the  brackets,  <  >,  denote  cross  sectional  average  values 

It  is  evident  that  for  adiabatic  flow  (a   =  0) 

w 

_  m  +  n  +  4 
o   m  +  n+2 


For  triangular  concentration  and  velocity  profiles 
(m  =  n  =  1) ,  C   reaches  a  maximum  of  1.5.   With  values  of  m 
and  n  between  2  and  3  parabolic  profiles  are  produced,  similar 
to  the  velocity  profiles  for  single-phase  laminar  flow. 
Higher  values  of  rri  and  n  (~7)  produce  profiles  similar  to 
turbulent  velocity  profiles.   As  a  flatter  profile  of 
concentration  and  velocity  is  reached,  C   approaches  a  value 
of  unity. 
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Zuber  and  Findlay  have  also  shown  that  the  weighted  mean 
gas  velocity  V   may  be  expressed  in  terms  of  C   and  a  weighted 
mean  gas  drift  velocity  V  . 


V  =  -43_ 

g    <a> 


<aV  > 

__S_  =  c  <j>  +  V  . 
<a>      o  J      gj 


where  V  .  =  V  -j  is  the  local  drift  velocity  of  the  gas, 

g-j  g     J  J.  J 


<j> 


Q     +  Qe 
_g__   f 

A 


<j  >  =  Q  /A 
Jg      g' 


and 


V 


93 


<aV  .> 
__21_ 

<a> 


Hence,  C   can  be  determined  from  the  slope  of  a  V   versus 
o  c  g 

<j>  plot,  and  the  weighted  mean  drift  velocity  from  the 
intercept  with  the  ordinate  axis.   A  typical  plot  in  the 


<J> 
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Hence  it  can  be  seen  that  plotting  V   versus  <j>  over  a 
wide  range  of  flow  conditions  will  yield,  at  least  for  round 
tubes,  some  information  on  the  void  and  velocity  profiles  in 
the  flow  duct.   Additional  information  on  concentration 
profiles  could  be  provided  by  a  traversing  void  fraction  probe, 
if  necessary. 

The  role  of  mixing  length,  I ,  in  two  phase  turbulent 
mixing  is  of  importance  if  one  considers  the  net  mixing  mass 

flow  rate,  w1  ,  ,  as  a  function  of  a  void  fraction  eddv 

net  ■* 

c\.  a 
diffusivity,  z    ,  and  a  concentration  gradient,  -5—. 

J    a     ,  3  dy 

_  net      1 1    1 1 

G__  .  __  .   .   =  — ; =  p.U..-p.U.. 


mixing 


[p-U-a.J+p  a.]  l~ 

f   3      g  d   dy 


[p_ (l-a. )+p  a . ]  £-g— 
Kf     1   *g  1    dy 


f   g    1   j    dy 


(p.e~P  )  (a- -a  . )  [£   3—] 

Kf  g   1  ]     ^y 

£ 

,      .da 

(p^-p  )  e  -5— 

f   g    ady 


where 


.  „2  du 
a      dy 


and 


u.  .  *=  appropriate  measure  of  the  turbulent  mean  flue- 
1  j  '     L  r 

tuating  velocities  of  liquid  and  gas,  here 
assumed  to  be  equal  (FT/HR) 
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G  .  .    =  net  mixing  mass  velocity  (LB  /HR  FT") 
mixing  ^  ■*     m' 

y  =  direction  perpendicular  to  the  main  flow 

between  adjacent  rods 

u  =  local  axial  velocity 

The  subchannel  densities  (p.,  p.)  here  refer  to  average 

representative  values  in  the  two  subchannels. 

With  knowledge  of  the  mixing  flow  rates  and  void  fraction 

profiles  it  may  be  possible  to  determine  at  least  the  order  of 

magnitude  of  I ,  the  two-phase  mixing  length.   Although  several 

investigators  have  proposed  that  G  .  .    is  independent  of 
^  r      £  mixing        ^ 

rod  gap  spacing,  c,  and  dependent  rather  on  hydraulic  diameter, 
other  investigators  have  proposed  the  opposite .   Lahey  and 
Schraub  (11)  discuss  this  matter  at  length.   The  order  of 
magnitude  of  the  mixing  length,  I,    may  shed  some  light  on  this 
controversy. 

While  the  above  formulation  seems  intuitively  correct, 
there  are  indications  that  it  is  probably  true  only  under 
certain  conditions.   One  such  condition  might  be  the  case  of 

two  identical  subchannels  in  communication.   In  this  case 

dcx 
when  the  void  fractions  are  equal,  (-=—  =  0)  ,  one  would  expect 

zero  net  mixing. 

A  broader  view  of  this  problem,  however,  suggests 

viewing  the  two  subchannels  in  communication  as  a  single 

channel  progressing  toward  a  condition  of  equilibrium 

concentration  and  velocity  distributions. 
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When  one  considers  a  subchannel  geometry  such  as  geometry 

(a)  in  Figure  1 ,  the  equilibrium  void  fraction  distribution 

dot 
probably  would  not  be  one  with  -=—  =  0.   This  supposition  is 

supported  by  work  discussed  by  Lahey  and  Schraub  which  shows 

that  there  is  a  distinct  variation  in  equilibrium  void  fraction 

from  the  corner  subchannel  (geometry  (a) ,  Figure  1)  and  the 

centermost  subchannels  of  the  rod  bundle.   Such  evidence 

supports  the  concept  of  a  bubble  drift  mixing  effect, 

previously  mentioned  in  Chapter  I. 

If  p.     and  p.     refer  to  average  representative 

1 ,eq      Kj ,  eq  *  r 

densities  in  the  two  subchannels  under  equilibrium  conditions, 
then,  using  the  same  notation  as  before,  one  can  show 


G  .  .    =  (p. -p.    )  u .  .  -  (p. -p.    )  u .  . 
mixing     Kj  Kj,eq  '  ij '    Ki  Ki,eq  '  ij 

S  (P  -p  )  E  (%2  .  da 

f   g    a   dy    dy 


Despite  the  fact  that  the  influence  of  void  fraction  on 
turbulent  mixing  is  not  well  understood,  it  would  seem  that  a 
first  step  in  any  such  investigation  is  the  determination  of 
average  void  fraction  over  a  wide  range  of  flow  conditions 
for  a  given  geometry. 

There  are  several  experimental  methods  for  measuring 
average  void  fraction.   The  simplest  method  is  the  use  of 
quick  closing  isolation  valves.   If  the  valves  can  be  shut 
simultaneously  and  sufficiently  fast  to  trap  a  representative 
sample  of  the  two-phase  mixture,  the  trapped  liquid  can  be 
metered  to  give  accurately  the  average  void  fraction.   The 
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application  of  this  method  to  a  non-circular  duct  necessitated 
the  maintenance  of  the  flow  duct  geometry  throughout  the 
valves  in  order  to  avoid  perturbation  of  the  flow. 
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CHAPTER  III 


Description  of  Apparatus 


General 


The  apparatus  utilized  to  determine  the  required  data 
was  a  basic  air-water  system  with  filters,  pressure  regulators 
and  surge  tanks  for  air  and  water,  mixing  devices,  test 
sections,  air-water  separators  and  flow  measuring  devices. 
Other  instrumentation  provided  were  air  and  water  temperature 
and  pressure  measurement,  pressure  drop  measurement,  and 
electric  flow  regime  detection.   The  void  fraction  determination 
was  accomplished  by  quick  closing  valves  and  read  directly 
through  a  transparent  test  section  of  lucite.   Figure  3  is  a 
schematic  representation  of  the  overall  experimental  apparatus. 
Letters  in  parentheses  throughout  this  section  represent _the 
various  portions  of  the  apparatus  as  indicated  in  Figure  3 . 

Figure  4  is  a  photograph  showing  the  major  portion  of  the 
apparatus.   On  a  support  framework  consisting  of  Dexion 
strengthened  by  plywood  was  mounted  the  three  foot  entrance 
section  (Q)  and  the  five  foot  interchangeable  test  sections  (S) , 
simulating  reactor  subchannels.   The  filters  (B,F) ,  pressure 
regulators  (C/G)  and  surge  tanks  (D,H)  were  located  on  the 
lower  level  with  the  mixing  devices  (N)  located  directly 
under  the  lucite  entrance  section. 

At  the  exit  of  the  vertical  test  sections  (V)  and  on  the 
upper  level  of  the  laboratory,  two  identical  air-water 
separators  were  mounted,  one  separator  for  each  subchannel. 
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Key  to  Figure  3 

A  -  Air  root  value 

B  -  Air  Filter 

C  -  Air  pressure  regulator 

D  -  Air  surge  tank 

E  -  Water  root  value 

F  -  Water  filter 

G  -  Water  pressure  regulator 

H  -  Water  surge  tank 

I  -  Air  throttle  valves 

J  -  Water  throttle  valves 

K  -  Water  inlet  rotameters 

L  ~  Air  inlet  thermocouples 

M  -  Water  inlet  thermocouples 

N  -  Mixing  devices 

0  -  Bypass  relief  valve  (in  use  for  void  fraction  testing  only) 

P  -  Bypass  valve  (in  use  for  void  fraction  testing  only) 

Q  -  Three  foot  entrance  section 

R  -  Electric  conductivity  (flow)  probes 

S  -  Five  foot  test  sections 

SA  -  Mixing  test  section 

SB  -  Void  fraction  test  section---right  channel 

SC  -  Void  fraction  test  section--lef t  channel 
T  -  Exit  blocks' 

TA  --  Mixing  section  exit  block 

TB  -  Void  fraction  section  exit  block 
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U  -  Pressure  taps 

V  -  Air-water  separators 

W  -  Water  spill  receiving  tank 
X  -  Outlet  air  flow  rotameters 

Y  -  Outlet  air  thermocouples 
Z  -  Outlet  water  rotameters 
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Photograph  of  Overall  Apparatus 
Figure  4 
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Air  flow  measurement  was  accomplished  by  use  of  rota- 
meters (X)  at  the  exit  of  the  air  from  the  separators.   The 
air  flow  meters  exhausted  to  the  atmosphere.   Because  of 
the  wide  range  of  air  flow  to  be  measured,  a  bank  of  flow- 
meters with  changeable  floats  was  provided  to  meet  the 
anticipated  flows. 

The  Martinelli-Nelson  correlation  for  void  fraction 
versus  quality  was  used  to  determine  the  steam-water  void 
fraction  corresponding  to  the  maximum  local  quality  existing 
in  a  boiling  water  reactor  at  1000  psig.   Using  this  same 
void  fraction  for  air-water,  the  Lockhart-Martinelli  (5)  and 
Wallis  (4)  correlations  were  used  to  predict  anticipated  air 

flow  rates  for  a  maximum  water  mass  velocity  of  1  x  10  LBm/ 

2 

HRFT  .   Because  of  the  large  number  of  flowmeters  required, 

air  flow  measurement  was  initially  installed  only  at  the  exit 
of  the  separators.  Thermocouples  (Y)  were  installed  prior  to 
the  air  flow  rotameters. 

Water  flow  measurement,  also  by  rotameters,  was  provided 
both  prior  to  the  mixing  devices  (K)  and  at  the  exit  of  the 
separators  (Z) . 

From  the  beginning  it  was  recognized  that  the  inaccuracy 
of  the  rotameters  (±.5%  repeatability)  might  present  problems 
in  determination  of  the  small  mixing  flows  but  would  be 
adequate  for  void  fraction  determination.   The  advantages  of 
a  continuous  rather  than  an  intermittent  device  and  the  high 
flow  rates  anticipated  for  the  higher  void  fractions  led  to 
the  decision  to  utilize  flow  meters  rather  than  an  intermittently 
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operated  receiving  tank  for  flow  measuring.   Data  using  the 
mixing  test  section  indicates  however,  that  an  intermittent 
flow  measuring  device  will  indeed  be  required,  primarily 
because  of  fluctuations  in  the  air  flow  at  very  low  and  very 
high  flow  rates. 

Separators 

The  separators  themselves  used  an  inverted  burette 
principle.   The  air-water  mixture  emanating  from  the  test 
sections  was  introduced  into  an  inverted  "can"  in  a  larger 
spill  tank  (see  Figure  5) .   The  level  in  the  outer  spill  tank 
was  kept  constant  by  the  spillway,  while  the  level  in  the 
inner  can  varied  with  air  flow  rate*   The  air-water  mixture 
was  separated  by  gravity  and  by  inertial  forces  by  use  of 
baffles  in  the  top  portion  of  the  inner  "can".   The  water 
spilled  into  a  spill  collection  tank  from  which  it  was  routed 
to  an  exit  water  flow  meter.   Slight  variations  in  the  water 
level  in  the  inner  can  caused  by  the  separation  process  were 
thought  to  be  the  cause  of  the  air  flow  fluctuations. 

Mixing  Device  and  Air  and  Water  Supplies 

The  mixing  device  (N)  is  shown  in  detail  in  Figure  6. 
The  air  was  introduced  radially  into  the  water  stream  by  a 
small  pipe  closed  at  the  end  and  drilled  with  small  holes 
along  its  length  and  circumference.   The  axis  of  the  pipe  was 
coincident  with  the  axis  of  the  water  pipe  and  mixed  flow  pipe. 

The  air  supply  chosen  was  the  200  psig  supply  in  the 
Engineering  Projects  Laboratory.   Since  the  mixing  loop  was 
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the  only  load  on  this  line,  the  supply  was  subject  to  a 
minimum  of  fluctuations.   It  was  found,  however,  that  because 
of  cycling  of  the  compressor,  it  was  necessary  to  bleed  air 
from  the  surge  tank  for  proper  control  by  the  air  pressure 
regulator  at  small  flow  rates. 

The  v/ater  supply  (street  water  pressure)  was  also  subject 
to  some  small  fluctuations  despite  use  of  the  regulator  and 
surge  tank.   Data  was  taken  largely  at  night  to  minimize  these 
fluctuations.   A  water  pump  may  be  installed  for  future 
operation  if  higher  pressures  and  less  fluctuations  are 
desired. 

Filtering  was  performed  with  a  Cuno  Micro-Kleen  Filter 
for  air  and  Fulflo  Honeycomb  Filter  for  water. 

Thermocouples  (L,'i)  were  installed  after  the  exits  of 
the  air  and  v/ater  throttle  valves  to  monitor  air  and  water 
temperatures  prior  to  mixing.   It  was  determined  that  these 
temperatures  would  best  represent  the  air  and  water  conditions 
in  the  test  section.   Only  slight  air  temperature  differences 
between  the  value  before  mixing  and  after  separation  verified 
this  assumption. 

Subchannel  Sections 

The  test  section  and  preceding  entrance  section  were 
machined  in  halves  from  lucite  using  a  special  milling  cutter. 
The  entrance  section  and  void  fraction  test  section  halves 
were  joined  with  ethelene  dichloride  (plexiglass  cement)  to 
produce  a  strong  joint  capable  of  withstanding  the  higher 
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pressures  anticipated  upon  the  closing  of  the  void  fraction 
isolation  valves.  The  mixing  section  halves  were  joined  by 
bolts  and  sealed  by  0  ring  material. 

The  entrance  section  and  interchangeable  test  sections 
were  joined  together  and  to  the  inlet  and  exit  blocks  with 
rubber  gaskets. 

The  three  foot  entrance  section  was  preceded  by  bypass 
connections  (0,P)  for  use  in  the  void  fraction  testing 
configuration.   For  least  disturbance  in  the  flow  pattern 
the  bypass  connections  were  placed  here  rather  than  just 
prior  to  the  isolation  valves.   Thus  a  three  foot  long  section 
identical  to  the  test  sections  in  geometry,  was  provided  for 
establishment  of  fully  developed  flow  prior  to  the  test 
sections  themselves . 

An  electric  flow  regime  detector  or  flow  probe  (R)  was 
installed  in  each  subchannel.   The  location  of  the  detector 
was  approximately  four  inches  before  the  outlet  of  the  entrance 
section.   With  the  location  of  the  flow  probe  in  the  entrance 
section/  flow  regime  detection  was  possible  for  all  test 
section  installations.   The  flow  probe  itself  was  manufactured 
using  a  24  gauge  stainless  steel  spring  wire  with  a  teflon 
shield.   A  silver  epoxy  paste  was  painted  on  the  lucite  sub- 
channel wall  for  the  wall  electrode.   Wire  leads  were  imbedded 
in  the  subchannel  wall  underneath  the  silver  epoxy  to  provide 
a  current  path.   A  detailed  description  of  the  electric  flow 
probe  installation  is  given  in  Appendix  B. 
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Entrance  Section 

The  entrance  section  was  constructed  as  indicated  in 
Figure  7.   A  1/32  inch  brass  strip  was  utilized  as  a  subchannel 
separator,  and  this  strip  was  sealed  in  a  machined  groove  in 
the  entrance  section  wall.   It  would  have  been  desirable  to 
use  a  thinner  strip  in  order  to  model  more  accurately  the 
conditions  in  the  actual  reactor  rod  bundle.   However,  two 
considerations  dictated  1/32  inch  as  the  minimum  strip  thick- 
ness.  Adequate  strength  under  the  possible  pressure  excursions 
from  valve  closing  in  the  void  fraction  test  section  required 
a  sufficiently  strong  strip.   Secondly,  the  sealing  of  the 
bottom  of  the  window  in  the  mixing  test  section  required  the 
machining  of  a  groove  in  the  fixed  (bottom)  portion  of  the 
strip  to  mate  with  a  knife  edge  in  the  upper,  movable  section 
of  the  strip.   It  was  this  upper,  movable  portion  of  the  strip 
which  allowed  the  adjustment  of  a  variable  sized  window 
between  the  two  subchannels  for  determination  of  turbulent 
mixing  and  for  complete  isolation  of  the  two  subchannels  from 
each  other  when  desired. 

Since  the  influence  of  hydraulic  diameter  on  mixing  was 
also  of  concern,  a  "plug"  was  provided  for  the  right  hand 
subchannel  in  the  entrance  and  mixing  test  sections  in  the 
form  of  a  brass  flat  bar  machined  to  specifications  and 
anchored  by  screws  to  the  subchannel  wall.   With  one  set  of 
brass  bars  in  place,  both  the  left  and  right  subchannels  were 
identical.   With  a  second  set  of  brass  bars  in  place,  the  right 
subchannel  was  significantly  larger  than  the  left  (see  Figure  7) 
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Mixing  Test  Section 

The  mixing  test  section  was  divided  into  upper  and  lower 
sections  of  three  feet  and  two  feet  respectively,  with  a 
total  length  of  five  feet.   The  bottom  of  the  window  was 
located  2  1/2  feet  from  the  bottom  of  the  mixing  test  section. 
The  upper  brass  strip  was  movable  to  allow  adjustment  of  the 
window  opening  and  was  withdrawn  through  a  special  packing 
gland  in  the  exit  block.   A  seal  on  the  movable  and  fixed 
brass  strips  was  maintained  by  C  clamp  pressure  on  the  exterior 
of  the  test  section.   Figure  8  shows  the  mixing  test  section. 

Void  Fraction  Test  Sections 

A  photograph  of  the  void  fraction  apparatus  is  shown  in 
Figure  9.   The  construction  of  the  void  fraction  and  entrance 
section  from  lucite  necessitated  a  rugged  mounting  arrangement 
for  the  void  fraction  isolation  valves.   The  void  fraction 
test  sections  apparatus  consisted  of  two  1  1/4  inch  ball  valves 
mounted  on  a  19  inch  x  60  inch  x  1/4  inch  thick  carbon  steel 
plate.   The  plate  in  turn  was  bolted  on  the  Dexion  framework. 
Valve  supports  were  provided  by  building  a  boxlike  framework 
of  1/2  inch  thick  steel  and  5/8  inch  thick  steel  bored  to 
accommodate  the  valve  bodies.   The  valve  bodies  were  machined 
at  the  inlet  and  outlet  for  brass  sleeves  which  were  soldered 
into  the  valve  inlets  and  outlets.   The  sleeves  served  two 
purposes:  1)  to  assist  in  fastening  the  valves  to  the  boxlike 
framework  and  2)  to  permit  insertion  into  the  valve  ends  of 
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cylindrical  lucite  valve  inserts  with  concave  ends  having 
the  same  subchannel  geometry  as  the  test  section.   A  similar 
cylindrical  insert  with  convex  ends  was  provided  for  each 
valve  ball.   Thus  the  subchannel  geometry  was  maintained 
throughout  the  entrance  section,  isolation  valves  and  void 
fraction  section.   The  only  interruption  in  the  subchannel 
geometry  were  the  -.020  inch  gaps  at  the  valve  balls  to 
preclude  hangup  of  the  ball  inserts  on  the  Buna  N  valve  seats. 

The  approximately  4  3  1/2  inch  isolation  section  between 
valve  balls  permitted  void  fraction  determination  from 
approximately  3  -  97%.   The  0-3  and  97  -  100%  range  of  a 
was  not  capable  of  being  read  directly  because  of  its  being 
hidden  by  the  valve  bodies. 

The  portion  of  the  test  section  between  the  valves  was 
provided  with  three  joints  sealed  by  O  rings.   Two  of  the 
joints  were  located  at  the  juncture  with  the  valve  inserts 
and  the  third  joint  was  located  in  the  center  of  the  section. 
Limitations  on  the  traverse  of  the  milling  machine  required 
the  test  section  between  the  isolation  valve  inserts  to  be 
divided  into  two,  approximately  20  inch  long/ sections . 

The  void  fraction  test  section  apparatus  consisted  of  a 
left  subchannel  and  a  right  subchannel  section,  each  capable 
of  being  installed  while  the  main  steel  mounting  plate  was  in 
place  on  the  Dexion  structure.   Thus,  there  were  interchangable 
valve  inserts,  valve  ball  inserts  and  void  fraction  isolation 
sections  for  both  the  left  and  right  hand  subchannels.   A 
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photograph  of  these  pieces  is  shown  in  Figure  10.   Figure  11 
shows  the  valve  ball  insert  in  place. 

Void  fraction  divisions  from  3  -  97.5%  were  scribed  on 
the  front  only  of  the  isolation  sections  by  use  of  a  milling 
machine.   Visual  error  in  reading  void  fractions  due  to 
refraction  was  found  to  be  minimal  since  the  meniscus  provided 
a  gauge  for  determining  when  the  line  of  sight  to  the  air- 
water  interface  was  above  or  below  the  line  of  sight 
perpendicular  to  the  interface. 

Pressure  taps  for  the  pressure  level  and  pressure  drop 
determination  were  installed  three  feet  apart  in  the  isolation 
section  of  the  void  fraction  test  sections.   Pressure  trans- 
ducer and  transducer  readouts  were  utilized.   This 
instrumentation  permitted  accurate  determination  of  the 
pressure  level  and  pressure  drop. 

The  valve  closure  device  selected  was  an  air  piston 
operated  mechanical  linkage  actuated  by  an  electric  solenoid 
simultaneously  with  a  similarly  operated  bypass  valve.   The 
main  isolation  valves  were  mechanically  linked  by  a  steel 
operating  rod  and  operated  by  a  single  air  piston  and  cylinder 
The  bypass  valve  was  operated  by  a  similar  linkage  but  with  a 
separate  air  piston  because  of  the  large  distance  from  the 
test  section  to  the  bypass  valve.   The  bypass  valve  (and 
operating  mechanism)  was  mounted  on  a  separate  Dexion  and  wood 
framework  on  the  lower  level. 

Each  air  operated  piston  and  cylinder  was  actucited  by 
separate  solenoids  connected  electrically  in  parallel. 
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Actuation  of  the  solenoids  closed  the  isolation  valves  and 

simultaneously  opened  the  bypass  valve.   Careful  adjustment 

of  air  operating  pressure  to  the  separate  cylinders  gave  a 

reliable  simultaneous  operation  in  70  milliseconds.   Preliminary 

calculations  using  the  Wallis  correlation  (4)  for  air  and 

water  showed  this  to  effectively  isolate  the  section  in  ~30% 

of  the  transit  time  for  the  gas  phase  at  80%  void  fraction 

and  G       =  1x10   lb/ft  hr .   The  channel  was  effectively 
water  J 

isolated  in  1/8  turn  of  the  isolation  valves  and  completely 
isolated  in  1/4  turn.   Valve  closing  times  were  monitored  on 
a  dual  beam  storage  oscilloscope  using  an  electrical  circuit 
with  potentiometers  mounted  on  the  valve  stems. 
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CHAPTER  IV 
Procedure 

The  right  subchannel  void  fraction  test  section  complete 
with  valves  and  inserts  was  mounted  first  above  the  right  hand 
subchannel  of  the  entrance  section.   By  so  doing,  the  left 
subchannel  of  the  entrance  section  was  blocked  off.   The  bypass 
valve  was  connected  to  the  right  subchannel  bypass  connection. 

The  inlet  v/ater  flow  rate  was  set  at  predetermined 
superficial  liquid  mass  velocities  matching  the  range  of  mass 

velocities  in  a  typical  boiling  water  reactor.   Six  superficial 

ft  ft  "? 

mass  velocities  ranging  from  0.5  x  10   to  1.0  x  10   LB  /HR  FT 

(with  a  0.1  x  10   increment)  were  investigated. 

For  each  liquid  mass  velocity,  the  air  flow  rate, 

monitored  at  the  exit  of  the  separators,  was  varied  throughout 

the  flow  range  of  the  air  flow  meter  bank.   Several  void 

fraction  readings  (at  least  three)  were  taken  at  each  air  flow 

rate  and  from  this,  an  average  void  fraction  was  determined. 

* 
It  was  found  that  in  the  bubbly  flov;  regime  there  was  little 

scatter  in  the  data,  while  in  the  slug  flow  regime  as  many  as 

five  or  six  trapped  samples  were  necessary  to  achieve  a  smooth 

void  fraction  curve.   Under  identical  conditions  individual 

void  fraction  readings  varied  by  as  much  as  8%  in  the  slug 

flow  regime. 

For  each  flow  condition  air  temperature  was  recorded  both 

prior  to  mixing  with  water  and  after  separation.   Water 

temperature  was  monitored  prior  to  mixing  only. 
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Both  pressure  level  and  pressure  drop  in  the  test  section 
were  recorded  for  each  point.   A  low-pass  filter  was  provided 
between  the  Ap  transducer  readout  and  the  vacuum  tube  volt- 
meter to  dampen  the  effect  of  the  pressure  drop  oscillations 
in  the  slug  flow  regime. 

The  electrical  conductivity  probe  and  visual  observations 
were  used  to  determine  the  flow  regime  boundaries.   Polaroid 
photographs  of  the  flow  probe  traces  were  taken  for  permanent 
record  of  the  flow  conditions.   Several  representative  photo- 
graphs are  shown  in  Appendix  B,  Figure  B-2.   These  photographs 
enable  the  determination  of  typical  slug  lengths. 

It  was  found  that  the  design  of  the  separators  was 
limiting  the  void  fraction  to  about  65%.   At  these  high  air 
flow  rates,  increased  pressure  drop  in  the  exit  piping  and  exit 
air  rotameters  caused  a  buildup  of  pressure  in  the  inner, 
inverted  "can"  of  the  separators.   This  pressure  exceeded  the 
hydrostatic  pressure  maintained  by  the  water  in  the  outer  drum 
and  resulted  in  escape  of  the  air  from  the  bottom  of  the  inner 
can.   Redesign  or  modification  of  the  separators  will  be 
necessary  to  reach  higher  void  fractions. 

When  the  required  data  had  been  obtained  for  the  right 
hand  channel,  the  corresponding  void  fraction  test  section 
was  removed  and  the  operations  were  repeated  for  the  left 
subchannel . 

Alteration  of  the  valve  closing  speed  over  a  range  of 
100-160  milliseconds  showed  no  change  in  the  average  void 
fraction  measured  for  typical  conditions  in  the  slug  flow 


regime.   Although  initial  valve  closure  times  were  of  the 
order  of  70  milliseconds,  an  average  of  100  milliseconds  vzas 
maintained  throughout  the  data  taking  period. 

The  instrumentation  was  calibrated  as  follows: 

a)  air  rotameters:   initial  calibration  with  known  air 
volume  flow  rates  using  a  gasometer 

b)  water  rotameters:   initial  calibration  with  known 
water  flow  rates  using  a  weigh  tank 

c)  pressure  level  transducer:   calibration  prior  to 
each  data  taking  run  using  a  mercury  manometer 

(normally  daily) 

d)  Ap  transducer:   electrical  calibration  prior  to  each 
data  taking  run  (normally  daily) 

e)  air  and  water  thermocouples:   no  calibration 
It  is  emphasized  that  throughout  the  range  of  flow 

conditions  covered  by  this  experiment,  no  attempt  was  made 
to  maintain  pressure  level  constant  at  any  location  in  the 
test  section. 
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CHAPTER  V 


Results 


The  experimental  results  of  this  thesis  are  presented 

graphically  in  this  chapter  for  the  highest  and  lowest  values 

of  superficial  liquid  mass  velocity  (1.0  x  10   and  0.5  x  10 

LB  /HR  FT2) . 
m 

Data  are  presented  in  tabular  form  for  these  superficial 
mass  velocities  and  intermediate  values  in  Appendix  C. 

Void  fraction  data  are  plotted  against  W  ,  the  air  mass 
flow  rate,  which  is  directly  proportional  to  quality  for  the 
low  air  mass  flow  rates  investigated.   Flow  regime  boundaries 
are  indicated  on  these  graphs.   Figures  12-15  show  these  data 
and  compare  them  to  correlations  by  Griffith  (12)  ,  Lockhart- 
Martinelli  (5)  and  Cravarolo  and  Hassid  (13).   The  Griffith 
prediction  applies  to  the  slug  flow  regime  only  and  is  plotted 
for  a  round  tube  having  the  same  hydraulic  diameter.   The 
Lockhart-Martinelli  correlation  was  developed  for  air  and 
water  and  many  other  liquids  in  round  tubes  and  horizontal 
adiabatic  flow.   Cravarolo  and  Hassid 's  data  was  obtained  for 
argon  and  water  or  ethyl  alcohol  in  several  geometries:   tubes, 
annuli,  and  a  seven  rod  bundle  in  a  circular  shroud. 

The  pressure  gradient  v/as  plotted  against  void  fraction, 
since  this  information  is  considered  to  be  of  interest  in 
establishing  the  proper  conditions  for  the  forthcoming  mixing 
experiments.   A  prediction  by  the  Lockhart-Martinelli  corre- 
lation for  total  pressure  gradient  versus  void  fraction,  is 
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included  for  comparison.   Measured  local  conditions  (i.e., 
pressure  level  and  temperature  at  the  level  corresponding  to 
the  window  location  in  the  mixing  test  section)  were  used  in 
the  Lockhart-Martinelli  correlation.   It  is  possible,  however, 
to  predict  void  fraction  and  total  pressure  drop  using  this 
correlation  without  knowledge  of  the  pressure  at  the  window 
location.   This  could  be  done  by  iterative  calculations  of 
pressure  gradient  and  void  fraction.  •  Pressure  gradient  data 
is  shown  in  graphical  form  in  Figures  16-19. 

Average  slip  ratio  is  shown  in  Figures  20-23.   The 
Marchaterre-Hoglund  correlation  (14)  for  slip  is  shown  in 
these  figures  for  comparison.   This  correlation  v/as  developed 
to  predict  velocity  ratios  in  vertical  channels  of  boiling 
water  reactors.   These  researchers  investigate  several  geo- 
metries in  formulating  their  correlation  (which  makes  no 
distinction  for  differing  geometry) .   Both  adiabatic  and 
diabatic  cases  were  studied  for  steam  and  water  in  rectangular 
and  circular  channels.   Air-water  experiments  were  conducted 
in  rectangular  channels  only.   Rod  bundle  geometries  were  not 

investigated. 

<j  > 
Average  gas  velocity,  V   =    *  ,    was  plotted  against 

average  volumetric  flux  density,  <j>,  in  Figures  24-27,  in 
the  form  set  forth  by  Zuber  and  Findlay  (9).   These  plots 
permitted  the  determination  of  C  ,  the  velocity  and  concen- 
tration distribution  parameter,  and  V  .,  the  weighted  average 

gas  drift  velocity.   In  all  cases  investigated,  V  .  was  found 

93 
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to  be  very  small,  less  than  1  FT/SEC.   Values  of  C   of  about 

1.25  were  obtained  for  a  superficial  liquid  mass  velocity, 

Gr:,    of  0.5  x  10   in  both  subchannels.   For  the  higher  mass 

velocity,  G,  =  1.0  x  10  ,  a  value  of  about  1.20  for  C   was 
J         f  o 

obtained. 

Intermediate  superficial  mass  velocities  for  the  left 
subchannel  will  be  investigated  in  -subsequent  work. 
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CHAPTER  VI 


Discussion  of  Results 


The  void  fraction  data  obtained  agree  quite  favorably 
with  the  various  predictions  plotted  in  Figures  12-15. 

The  Griffith  slug  flew  correlation  (12)  was  developed 
using  round  tubes,  annuli,  and  circular  shroud  bundle  geometries 

The  void  fraction  under  adiabatic  conditions  is  given  by: 

R   =  a  •■■■  -% 


g        Q  +Q.-J-V,  A 
g   f   b 


where 


Vb  =  ^/gD  +  V-^I) 


D  is  defined  differently  for  round  tubes  and  other  geometries, 
and  K,  and  K„  are  also  geometry  dependent. 

At  first  glance  it  would  appear  proper  to  utilize  Griffith's 
rod  bundle  curves  for  estimation  of  the  constants  K,  and  K_  in 
his  correlation.   Griffith's  work  entailed  the  use  of  a  2.0" 
shroud  diameter  bundle  with  seven  rods.   Keeping  the  shroud 
diameter  constant,  he  varied  the  rod  diameter,  thereby  altering 
the  hydraulic  diameter.   The  bubble  velocities  Griffith  ob- 
served were  bundle  averages,  not  the  average  subchannel 
velocities  of  interest  in  the  subchannels  tested  for  this 
thesis.   Hence  the  asymptotic  values  of  K,  and  K„  for  round 
tubes  were  used  in  applying  Griffith's  correlation.   When  used 
in  this  manner,  Griffith's  correlation  gave  results  which 
agreed  very  favorably  with  the  data. 
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For  the  right  subchannel,  Griffith's  correlation  for 
round  tubes  produced  void  fractions  somewhat  less  than  those 
measured.   In  the  left,  more  restricted  subchannel,  the 
Griffith  prediction  yielded  a  higher  void  fraction  than  that 
observed.   Intuitively,  one  would  anticipate  generally  lower 
void  fractions  here  than  in  round  tubes  because  of  retardation 
of  the  liquid  in  sharp  corners,  etc.   The  fact  that  this  same 
effect  was  not  observed  in  the  right  subchannel  is  subject  to 
speculation. 

Comparison  of  Griffith's  slug  flow  correlation  with 
Zuber  and  Finley's  (9)  formulation  gives  the  following 
relationships : 

«  = ?* • 


Qf+Q„ 

Q  +Qf+AK1/5D+K2  <-a-*)A 


Qg/A 


Q.+Q 
(1+K2)  (-^-^-H^/gD 

or  in  Zuber  and  Findlay ' s  terminology 

<j  > 

<a> 


( l+K2"f^J>+K1/gb 


c  <j>+v  . 

o     gj 


hence 


C   =  1+Kn  and  V  .  =  K,  ZgD 
o       2       gj     1  3 


The  values  of  C   obtained  from  the  V   versus  <j>  plots  (Figures 

o  g  c  j 

24-27)  yield  values  of  K   of  approximately  0.20  and  0.25.   The 
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intercept  on  the  ordinate  axis,  V  ...  ,  is  difficult  to  specify 
exactly  but  seems  to  range  from  zero  to  about  0 . 5  FT/SEC . 

Use  of  the  asymptotic  values  for  circular  pipes  for  K, 
and  K„  and  equivalent  pipe  diameter  for  D  in  Griffith's 
correlation  agrees  very  closely  with  the  information  gleaned 
from  the  V  versus  <j>  plots,  corroborating  the  assumption 
that  Griffith's  correlation  for  cirular  pipes  applies  to  the 
geometry  investigated. 

Griffith's  correlation  used  in  this  manner  predicts 

values  for  C   of  1.20  and  V  .  of  0.3  4  and  0.39  FT/SEC  for  the 
o  gj 

left  and  right  subchannels  respectively. 

The  values  of  C   of  about  1.20  and  1.25  obtained  for  the 

o 

r 

two  liquid  superficial  mass  velocities  of  1.0  x  10   and 

0.5  x  10   LB  /HR  FT'  respectively  are  in  consonance  with  a 

flatter  profile  expected  from  doubling  the  water  mass  velocities 

If  Zuber  and  Findlay's  definition  of  C   for  round  tubes  is  used 

J  o 

and  the  concentration  and  velocity  profiles  are  assumed  to  be 

identical,  exponent  values  of  m  and  n  of  3(Gf=0.5  x  10   LB  /HR 

FT  )  and  4(G  =1.0  x  106  LB  /HR  FT2)  are  obtained  for  both 
£  m 

subchannels.. 

The  Lockhart-Martinelli  void  fraction  predictions  fall 
below  the  observed  data  in  the  region  of  intermediate  void 
fractions.   In  the  low  and  high  void  fraction  regions  the 
correlation  seems  to  agree  quite  well  with  the  data.   This 
effect  does  not  seem  to  be  flow  regime  dependent'  at  the 
lower  end  of  the  region  of  discrepancy,  since  it  is  present 
in  the  bubbly  flow  regimes  for  the  right  subchannel  at 
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G,  =  1.0  x  10   LB  /HR  FT"  and  in  the  slug  flow  regime  for  the 

f  m  3         J 

other  flow  conditions  analyzed.   As  the  annular  flow  regime  is 
approached,  the  Lockhart-Martinelli  prediction  agrees  more 
closely  with  the  data  and  may  even  give  higher  void  fractions. 
Slug-annular  flow  was  observed  in  the.  left  subchannel  at  the 
higher  gas  flow  rates  investigated,  but  complete  annular  flow 
was  never  reached.   The  fact  that  the  Lockhart--Martinelli 
predictions  fall  below  the  observed  void  fraction  data,  rather 
than  above  as  one  might  expect,  is  unexplained. 

The  Cravarolo-Hassid  correlation  (13)  for  steam-water 
is  included  in  the  void  fraction  versus  W   curves  for  a 

g 

comparison  with  their  correlation  and  for  steam-water  and  air- 
water  mixtures. 

The  pressure  gradient  versus  void  fraction  curves  (Figures 
16-19)  agree  quite  favorably  in  shape  with  the  curves  predicted 
by  the  Lockhart-Martinelli  correlation,  but  are  consistently 
below  the  predicted  values. 

Scatter  in  pressure  gradient  data,  particularly  at 
higher  void  fractions  in  the  slug  flow  regime  can  be  attributed 
to  two  sources  of  error: 

1)  pressure  drop  fluctuations  present  in  slug  flow 
conditions,  the  effect  of  which  could  not  be 
completely  eliminated  by  the  low  pass  filter,  and 

2)  fluctuations  in  liquid  mass  flow  rate  caused  by  poor 
control  by  the  water  pressure  regulator.   These, 
though  small,  have  a  significantly  higher  effect  on 
pressure  gradient  than  on  void  fraction  (such 
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fluctuations  were  more  prevalent  in  the  slug  flow 
regime) . 

Average  slip  ratio,  plotted  versus  gas  mass  flow  rate, 
W  ,  in  Figures  20-23,  seems  to  be  in  rough  agreement  with  that 
predicted  by  the  Marchaterre-Hoglund  correlation.   For  the 
right  subchannel,  the  observed  slip  ratios  fell  generally 
below  the  predicted  values.   On  the  other  hand,  observed  slip 
ratios  in  the  left  subchannel  were  consistently  greater  than 
predicted  values. 

Since  the  slip  ratio  is  related  to  void  fraction, 

s  =  v  /v  =   x  ■  a^i  m 

S    VVf    (l-X)     a    pg 

one  would  anticipate  generally  higher  slip  ratios  for  lower 
void  fractions  with  equal  quality  flows  in  different 
geometries . 

The  observed  slip  ratio  for  the  left,  more  restricted 
subchannel,  then  would  probably  be  higher  than  that  predicted 
for  a  more  open  (circular  or  rectangular)  geometry,  as  is  the 
case.   The  fact  that  the  data  for  the  right  subchannel  falls 
below  the  predicted  curve,  at  least  for  all  but  the  higher 
void  fraction  region,  seems  reasonable,  since  it  may  be  even 
less  restrictive  than  the  average  of  Marchaterre  and  Hoglund's 
geometries . 

Scatter  in  the  slip  ratio  data  may  be  attributed  to  the 
fact  that  the  accumulation  of  errors  in  the  void  fraction 
measurement  and  flow  measurement  become  magnified  in  this 
display  of  data.   This  is  particularly  true  at  lower  void 
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fractions,  where  a  small  error  in  void  fraction  is  multiplied 

beyond  proportion  in  the  slip  ratio.   The  fact  that  the  scatter 

is  not  uniform  about  a  slip  ratio  of  approximately  unity  at 

lower  void  fractions  may  be  indicative  of  inadequate  operation 

of  the  mixing  device  at  a  low  gas  mass  flow  rates. 

An  analysis  of  several  representative  electrical 

conductivity  probe  traces  in  the  slug  flow  regime  yielded 

bubble  lengths  varying  from  1/2  -  6  inches  with  as  much  as 

4  inches  between  bubbles.   This  data  was  obtained  for  the 

right  subchannel  at  G,  =  0.5  x  10   LB  /HR  FT   with  a  =  .66. 
3  f  m 

In  flow  conditions  such  as  these,  a  wide  scatter  in  void 
fraction  data  occurred.   The  test  section  length,  when  compared 
with  bubble  lengths,  was  probably  too  small  to  achieve  an 
accurate  statistical  average  in  a  single  sample. 

A  change  in  pressure  at  the  level  corresponding  to  the 
window  location  in  the  mixing  test  section  will  obviously 
affect  void  fraction.   An  experiment  to  determine  the  magnitude 
of  this  variation  was  conducted  by  approximately  doubling  the 
gage  pressure  at  the  window  location  for  flow  conditions  in 
the  bubbly  flow  regime.   A  decrease  in  void  fraction  of 
approximately  .04  was  observed. 

It  will  be  necessary  to  throttle  the  exit  of  either  or 
both  subchannels  in  the  mixing  experiments  to  produce  equal 
pressure  level  on  both  sides  of  the  window.  It  is  expected 
that  the  effect  on  void  fraction  caused  by  such  throttling  can 
be  predicted  using  the  Lockhart-Martinelli  correlation.  It 
will  be  necessary,  however,  to  substantiate  this  by  experimental 

results . 
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CHAPTER  VII 


Conclusions 


An  experimental  apparatus  was  designed,  constructed  and 
yielded  satisfactory  data  for  void  fraction  over  a  wide  range 
of  flow  conditions.   Higher  void  fraction  regions  beyond  those 

covered  (a     =  .70)  may  be  reached  with  modification  or 

max  ■* 

redesign  of  the  air-water  separators. 

The  data  obtained  agreed  quite  favorably  with  various 
correlations,  although  some  discrepancy  was  present.   Notably, 
the  void  fraction  data  were  generally  somewhat  higher  than  the 
Lockhart-Martinelli  predictions  in  the  intermediate  void 
fraction  range.   The  Griffith  slug  flow  correlation  agreed 
quite  well  with  the  data  in  most  cases. 

Pressure  gradient  varied  with  void  fraction  in  such  a 
manner  that  it  should  be  possible  to  match  pressure  level  and 
pressure  gradient  under  flow  conditions  of  widely  different 
void  fractions  in  the  mixing  experiment. 

Plots  in  the  V   versus  <j>  plane  using  Zuber  and  Findlay's 
formulation  gave  a  rough  indication  of  the  peaking  of  the 
velocity  and  concentration  profiles  for  varying  superficial 
liquid  mass  velocities.   Detailed  profiles  should  be  provided 
by  other  means,  if  needed. 

The  electrical  conductivity  probe  was  very  useful  in  the 
identification  .of  flow  regimes  and  provided  a  means  for 
estimating- slug  and  bubble  lengths.   The  design  of  the  wall 
electrode  for  a  normally  non-conducting  channel  wall  proved 
to  be  satisfactory. 
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Flow  conditions  in  the  subchannel  geometries  investigated 
are  probably  not  representative  of  rod  bundle  average  conditions , 
but  rather  representative  of  individual  subchannels.   Indeed, 
the  purpose  of  these  experiments  was  the  simulation  of  sub- 
channels.  Preliminary  indications  suggest  that  the  right 
subchannel,  which  models  the  flow  region  between  four  fuel 
rods,  is  representative  of  the  actual  geometry. 
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APPENDIX  A 
Two-Phase  Conservation  Equations  and  Void  Fraction 

A  better  understanding  of  the  influence  of  void  fraction 
on  tv70-phase  turbulent  mixing  may  be  achieved  through  examina- 
tion of  the  relevant  conservation  equations. 

Several  multichannel  computer  programs  have  been  developed 
using  these  equations  to  predict  flow  conditions  and  ultimately 
the  burnout  point  for  nuclear  power  reactors  with  ventilated 
parallel  channels.   The  COBRA  Code,  (7),  developed  for  single 
and  two-phase  flow,  was  selected  as  an  example  of  present 
attempts  to  solve  the  problem  of  mixing  as  applied  to  burnout 
prediction. 

The  assumptions  of  the  COBRA  Code  are  as  f ollows : 

a)  One-dimensional,  two-phase  flow  is  described  by  two 
average  velocities,  a  gas  velocity,  V  ,  and  a  liquid 
velocity,  V ,- .   The  velocity  or  slip  ratio,  S=V  /Vf, 
is  generally  larger  than  one. 

b)  A  static  pressure  drop  can  exist  between  adjacent 
subchannels,  resulting  in  diversion  cross  flow 

c)  In  diversion  cross  flow,  the  cross  flow  enthalpy  is 
that  of  the  subchannel  from  which  it  originates 

d)  Turbulent  mixing  can  transfer  enthalpy  but  not  mass 
Also  implicit  in  the  assumptions  is  the  representation  of 
subchannel  variables,  such  as  enthalpy,  velocity,  etc,  by 
single  average  values. 
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Lahey  and  Schraub  (11)  suggested  that  a  no  net  mass 
transfer  approach  cannot  describe  two-phase  turbulent  mixing 
adequately,  and  they  propose  an  equal  volume  interchange  as 
a  more  reasonable  assumption  from  mixing  length  theory.   This 
assumption  requires  experimental  verification.   The  small 
amount  of  mixing  data  collected  to  date  using  the  mixing 
apparatus  indicate  that  for  parallel  channels  subjected  to 
equal  pressure  gradient  there  may  be  'a  net  turbulent  mixing 
flov7  by  volume  and  by  weight. 

Another  point  raised  by  Lahey,  Shiralkar  and  Radcliffe  (3) 
is  that  the  cross  flow  enthalpy  is  not  alv/ays  the  same  as  that 
of  the  donor  subchannel.   One  explanation  for  this  condition 
may  be  a  different  velocity  for  each  of  the  two  phases  in  the 
diversion  cross  flow. 

An  example  with  two  parallel  and  communicating  subchannels 
has  been  selected  to  illustrate  the  use  of  the  conservation 
equations,  taking  into  consideration  the  effect  of  these 
proposed  modifications  to  the  equations  presented  in  reference 
(7). 

The  following  additional  assumptions  have  been  included 
in  these  equations : 

a)  A  net  turbulent  mixing  mass  flow  rate  can  exist  in 
two-phase  flow 

b)  The  diversion  cross  flow  liquid  and  vapor  velocities 
may  differ 

c)  An  axial  velocity,  u,  is  defined  in  order  to  conserve 
the  momentum  of  the  mixture: 
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U  =  V. (1-X)  +  V  X  =  Gv' 
f  g 


where  v'  is  defined  such  as  to  preserve  the  momentum  of  the 
mixture: 

2       y2 
f      hg 


v 


(1-X) 
(l-a)p 


The  following  mass  continuity,  energy  and  momentum 
equations  apply  to  the  control  volume  shown  below  for  the 
momentum  equation. 
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MASS  CONTINUITY  EQUATION 


dm. 

-= —  =  -w .  .  +  w  .  .  -  w .  .  '  +  w  .  .  ' 
dx       13     31     ij      -j± 


where  either  w. . t    w.  .  ,  or  both  are  zero  depending  upon  the 
1 J    D1 

existence  and  direction  of  diversion  cross  flow. 


w.  .  '  =  cG.  .  *  =  c  [  (1-a.  )  p..  u.  .  -  +a .  p   u.  .    ] 
13       ID  1  f  i]^1,  1  g1  ID  /9 

where  [u._.  _  |  and  |u.  .   |  are  appropriate  measures  of  the 
13  / 1         id  / g 

turbulent  components  of  the  liquid  and  gas  velocities 
perpendicular  to  the  main  axial  flow  (x  direction),  e.g., 
RMS  values  or  averages  of  the  absolute  values.   Changes  of 
the  void  fraction,  a.  ,  within  the  length  dx  are  neglected  in 
the  above  equation.   Similarly, 


w..1  =  cG . . '  =  c[ (l-a.)pr  u. .  _  +a . p   u.. 
D1       Di  3   f   Di/f   D  g1  D1! 


but   u .  .  rr=   u  .  .  ^  and   u .  .    =   u .  . 

1  13 ,f '   1  31, f 1     1  id  ,g'   '  31, g 


hence 


JJL 3_  =  q  ,  .   =  (a.  -a  .)  p^-lu.  .  r.  I  +  (a  .--a  .  )  p  lu  .  . 

c        mixing    1   ]  Kf  i],f    3   1   9"'  13  ,9 


The  diversion  cross  flow  terms,  allowing  for  slip  in 
the  direction  perpendicular  to  the  main  flow  become 


w .  .  =  c  [  ( 1  -  a  .  )  p  -u .  .  j-   +  a.p  u..   ] 
ij  1   r  ij,f     1  g  ID r9 


for  flow  from  subchannel  i  to  subchannel  j  and 
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w  .  .  =  c[(l-a.)p  u  .    +  a.  pi,..   ] 

for  flow  from  subchannel  j  to  i. 

The  appearance  of  the  void  i taction  difference  between 
adjacent  subchannels  in  the  turbulent  mixing  mass  flow  rates 
is  a  strong  indication  that  void  fraction  can  be  used  as  a 
correlating  parameter  for  mixing.   if  one  considers  the 
anticipated  influence  of  void  fraction  on  the  fluctuating 


mixing  velocities,   u  .     and   u , .   I   the  assumption  that 
void  fraction  is  an  important  pai ameter  in  mixing  seems  even 
more  conclusive. 

Such  a  correlation  need  not  necessarily  be  linear,  since 
(as  noted  by  Lahey  and  Schraub  (1  I ) )  turbulent  mixing  seems  to 
be  strongly  dependent  on  the  two-phase  flow  regime,  peaking 
somewhere  in  the  slug  flow  regime.   Judging  from  the  disturbance: 
caused  by  the  passage  of  bubbles  of  various  sizes,  this 
suggestion  seems  reasonable. 

ENERGY  EQUATION 

d  (m .  h  .  ) 

— d*     =  q,  -  w.  .h.  -I-  w  'h   --  w.  .  'h.  '  +  w.  .  'h.  ' 

Where,  as  before,  either  one  or  both  of  the  diversion  cross 
flow  terms  may  be  zero,  depending  on  the  existence  and  direction 
of  the  diversion  cross  flows. 

Assuming  a"  uniform  void  fraction  distribution,  the  mixing 
energy  transfer  terms  are 
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w..'h.'    =    c  [  (1-a  .  )  p-.h_  lu.  .     ,|+a.o   h    lu..       I] 

13     1  1  K f  f   1 3  , f '     i  g  g '   1 3  / g ' 


w  . 


..'h.'    =   c  [  (1-a  .  )  Oj-h-   u.  .     r.  M-a  .  p   h      u.. 

did  d     f  f '   iD/f      d  g  g1   iD.g1 


and  the  diversion  cross  flow  energy  transfer  terms  are 


w.  .  h  .  =  c  [  (1-a  .  )  p_h  _u  .  .  ,_  +  a.  p  h  u.  .   ] 
ID  i  i   f  f  ID rf  i  g  g  ID i9 


w..h.  =  c[.(l-  a. )  pJi_u  .  .  ,.  +  a  .  p  h  u  •  ■   ] 

DID  ]   f  f  ]i<f     D  g  g  31 ,g 


MOMENTUM  EQUATION 

gp 


g  [A.p.-F.dx-A.  (p.+dp.  ) A.dx] 

^C   1^1   1     l  L  l  l  l   q     l 


3 

c 


f  ^w  .  .  u  .  dx  +  f  mw.  .'  '  u  .  dx  +  m  .  u  .  +  d  (m  .  u  .  ) 
D  13  l       T  id   1       11       11 


(f^w..u.dx  +  f_w..'u.dx  +  m.u.) 
D  ]l  ]        T  31   3        11 


Here,  as  in  the  COBRA  Code,  f  and  f  account  for  non- 
uniform velocity  and  void  fraction  (concentration)  profiles 
in  a  subchannel.   The  static  average  density  is  defined  as 


p .  =  (1-a. ) p,  +  a . p 
Ki        if     i"g 


Expanding    the    axial  momentum   change    term 


d  (m .  u  .  )    =  m .  du  .  +u  .  dm  .  =m .  d  (G  .  v  .'  )  -i-u  .  (w  .  .  -w  .  .  +w  .  .  '  -w .  .  '  ) 
11  1      1      1      1      111         131      13       31         13 

=   m .  [ v ! dG . +G . dv ! ]    +   u .  (w . . -w .  . +w  .  .  ' -w .  .  ' ) 
11      1      1      1  131      13      31         13 
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2 

u  .  dm .  -f  A  .  G  /dv  !  +  u  .  (w  .  .  ~w .  .  -f-w  .  .  '  -w  .  .  '  ) 
1   1   13.   1     i   31   lj   ji  ij 


2 

A . G . dv !  +  2u . (w . . -w . . +w . . ' -w . . ' ) 
li   l      l   ji  ±j       ji    ij 


2  2 

F.    f.<j>„GT 
l     l Y  £  i 


A.    2g  p_D 

1  JcKf  e,i 

where  f.  is  the  friction  factor  obtained  assuming  that  all  the 
flow  is  liquid  and 

2  dp/dx)FR  Tp 
♦t  -  dp/dx)FR  £ 

is  the  two-phase  frictional  pressure  drop  multiplier. 
Hence  the  momentum  equation  can  be  reduced  to 

dp.    f.cj)2         0  G2    dv! 

dx     2g  prD   .    i  q     dx 

^cKf  e,i  ^c 

friction  pressure  drop      acceleration  pressure 

drop  due  to  boiling 

+p.  2~  + — I — [u.  (w.  .  (f  -2)+w.  .  (2u.-f^u.) 

i  g  g  A.    iv  13   D      31   .1   D  j 

gravity  pressure  momentum  exchange  terms 

drop 

+  u.w. , ' (f_-2)+w. . ' (2u.-f_u.) ] 
1  1 3    T      D1     1   T  3 

It  is  apparent  that  a  thorough  study  of  the  turbulent 
mixing  phenomenon  as  it  affects  energy  and  momentum  transfer 
must  consider  the  concentration  and  velocity  profiles.   The 
lumped  parameter  approach  may  be  justified,  but  the  actual 
determination  of  a  valid  mixing  length  will  be  dependent  upon 
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not  only  the  mean  void  fraction  and  mean  velocity  in  each 
channel  but  also  on  concentration  and  velocity  distributions 
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APPENDIX  B 
Descri p tion  of  the  Electrical  Conductivity  Probes 

The  electrical  conductivity  probes  were  mounted  near  the 
discharge  end  of  the  entrance  section  for  both  left  and  right 
subchannels.   In  this  location  the  probes  allowed  detection  of 
flow  regimes  for  each  of  the  three  test  section  configurations. 
Figure  B--1  shows  the  flow  probe  installation. 

The  most  difficult  problem  encountered  in  the  construction 
of  the  probe  was  to  provide  a  conducting  electrode  at  the  wall, 
since  the  wall  material  was  lucite  and  non-conducting.   Several 
methods  were  investigated,  and  a  conducting  silver  epoxy  paint 
was  selected.   This  paste-like  material  (EPO-TEK  418H)  was 
painted  in  a  very  thin  coat  on  the  lucite  walls  of  the  sub- 
channels for  a  total  subchannel  length  of  six  inches ,  centered 
at  the  flow  probe.   A  portion  of  the  conducting  wall  was  formed 
by  the  brass  strip  separating  the  two  subchannels.   In  the  right 
subchannel  an  additional  portion  of  the  conducting  v/all  was 
formed  by  the  "plug"  (brass  flat  bar)  inserted  in  order  to 
alter  the  hydraulic  diameter  of  this  subchannel.   Outside 
connection  with  the  wall  electrode  was  accomplished  by 
threading  a  small  mu It i strand  wire  through  a  hole  in  the 
lucite  wall  and  laying  several  strands  in  a  shallow  groove 
machined  in  the  subchannel  wall.   The  groove  was  then  filled 
with  the  silver  epoxy,  covering  the  wire  strands.   This 
produced  a  sufficiently  smooth  surface.   The  epoxy  was  cured 
by  careful  use  of  a  heat  lamp,  insuring  that  the  temperature 
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was  not  sufficiently  high  to  cause  plastic  flow  or  warping 
of  the  lucite. 

The  probes  themselves  were  cut  from  2  4  gauge  stainless 
steel  spring  wire,  rounded  at  the  probe  tip  and  shielded  with 
thin  wall  teflon  spaghetti,  except  for  an  approximately  1/16" 
bare  portion  at  the  tip.   The  teflon  shield  was  not  wet  by  the 
liquid  in  the  subchannel  and  consequently  provided  insulation 
between  the  probe  tip  and  the  wall  electrode.   The  overall 
probe  diameter,  including  the  teflon . shield ,  was  about  0.035", 
approximately  8%  and  10%  of  the  hydraulic  diameters  of  the 
right  and  left  subchannels  respectively.   While  a  smaller 
probe  would  have  been  more  desirable  insofar  as  perturbation 
of  the  flow  is  concerned,  it  would  not  have  had  sufficient 
resistance  to  bending  when  buffeted  by  the  flowing  air-water 
mixture.   A  packing  gland  (Conax  MIC-062A)  was  screwed  into 
an  external  mounting  block  and  provided  a  seal  for  the  probe 
as  it  passed  through  the  subchannel  wall. 

The  electrical  circuit  was  identical  in  principle  to 
that  described  in  reference  (15) .   A  six  volt  battery  was 
used  to  supply  power.   Shielded  cable  was  used  for  connections 
to  eliminate  noise. 

Typical  probe  traces  for  various  flow  conditions  are 
shown  in  Figure  B~2.   These  Polaroid  photographs  of  oscillo- 
scope traces  were  taken  for  each  of  the  data  points  recorded 
for  the  right  subchannel. 

In  explanation  of  the  photographs,  signals  near  the  top 
of  the  photographs  indicate  a  current  path  through  the  water. 
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Signals  at  the  bottom  of  the  photographs  indicate  interruption 
of  the  current  path  by  the  passage  of  an  air  bubble. 
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Right  Subchannel  G.  =  0.6  x  10   LB  /HR  FT 
J  f  m 


W      P  at    p   at   Flow   <a>   No.  of    dp   T         T 

window   ?  ,    Reqime      Samples  d.1        g  inlet   f  inlet 
window    ^  L 

loc.    , 

loc. 

LBm      .  _     LB„  psi     0T.       OT^ 


psia       m 


HR     ^         — f  FT 

FT 


15.6  29.92  0.1470  S  .673  7  .417  73.0  63.0 

10.78  27.23  0.1310  S  .658  10  .420  80.0  59.5 

8.45  26.20  0.1252  S  .631  13  .403  80.0  59.5 

6.18  24.99  0.1195  S  .628  8  .394  80.0  59.5 

4.11  23.15  0.1110  S  .559  6  .367  80.0  59.5 

2.350  21.41  0.1010  B-S  .518  3  .406  80.0  59.5 

2.03  21.64  0.1045  B  .500  3  .357  80.0    .  39.5 

1.-485  20.72  0.0990  B  .388  3  .350  80.0  59.5 

1.320  20.42  0.0979  B  .366  3  .350  80.0  59.5 

0.819  19.54  0.0917  B  .259  3  .360  80.0  59.5 

0.394  18.88  0.0885  B  .138  3  .387  80.0  59.5 

0.314  18.78  0.0880  B  .126  3  .387  80.0  59.5 

0.0996  18.50  0.0865  B  .065  3  .400  80.0  59.5 
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Right  Subchannel  G  =  0.7  x  106  LB  /HR  FT2 

W      P  at    pg  at   Flow   <a>   No.  of    dp   T   .  _  .  T_  .  .  , 

a  ■    .  •,     ■  -,     ~  ~        ^  «    9  inlet   f  inlet 

3    window  window   Regime      Samples   dx 

loc.    loc. 


LB 

m 

HR 

psia 

LB 
m 

FT3 

psi 

FT 

°F 

op 

15.72 

29.53 

0.1430 

S 

0.651 

6 

.453 

73.5 

63.5 

10.52 

27.76 

0.13  4  5 

S 

0.611 

6 

.4  25 

81.5 

61 

9.46 

28.62 

0.1375 

S 

0.612 

6 

.430 

81.5 

61 

7.92 

27.81 

0.1345 

S 

0.616 

13 

.423 

81.5 

61 

6.15 

28.82 

0.1390 

S 

0.586 

7 

.433 

81.5 

61 

4.12 

27.93 

0.1340 

B-S 

0.524 

■  9 

.457 

81.5 

61 

3.32 

24.41 

0.1170 

B 

0.522 

5 

.416 

81.5 

61 

1.81 

22.59 

0.1082 

B 

0.434 

3 

.384 

81.5 

61 

1.42 

20.79 

0.1045 

B 

0.342 

3 

.374 

81.5 

61 

1.00 

21.79 

0.0995 

B 

0.272 

3 

.377 

81.5 

61 

0.558 

'  19.80 

0.0950 

B 

0.160 

3 

.390 

81.5 

61 

0.252 

19.50 

0.0935 

B 

0.057 

3 

.403 

81.5 

61 
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Right  Subchannel  G  -  0.8  x  106  LB  /HR  FT2 

W     P  at    pg  at   '   "     <a>   No.  of  dp   T   .      Tr  . 

a     ...     ?  ,    „Flov/  n  %-tt   g  inlet   f  inlet 

^    window  window  Regime       Samples  d£ 

loc .    loc . 


LB 
m 

HR 

psia 

LB 

m 

FT3 

psi 

FT 

op 

°F 

14.50 

30.44 

•  14  7  5 

S 

0.628 

6 

0.493 

7  4.5 

63.5 

10.50 

29.58 

•  1405 

S 

0.609 

7 

0.463 

84.0 

59.5 

8.4  9 

28.98 

•1375 

S 

0.580 

7 

0.465 

84.0 

59.5 

5.73 

27.94 

.1340 

S 

0.541 

8 

0.493 

81.0 

59.5 

5.18 

27.22 

.1307 

s 

0.534 

7 

0.476 

81.0 

59.5 

4.33 

27.05 

.1292 

B-S 

0.527 

3 

0.433 

81.0 

59.5 

3.09 

25.52 

.1223 

B 

0.50  6 

3 

0.410 

81.0 

59.5 

1.852 

23.78 

.1140 

B 

0.396 

3 

0.390 

81.0 

59.5 

1.39  0 

22.38 

.1065 

B 

0.304 

3 

0.386 

81.0 

59.5 

1.052 

21.87 

.1042 

B 

0.248 

3 

0.381 

81.0 

59.5 

0.627 

20.98 

.1001 

B 

0.164 

3 

0.390 

81.0 

59.5 

0.242 

20.00 

.0955 

B 

0.052 

3 

0.400 

81.0 

59.5 
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Right    Subchannel      G    =    0.9    x    10      LB    /HR   FT 
3  m 


W  P    at         Pq    aC    Flov/         <a>      No.    of         dp      T       .     .     .     T^    .     _     . 

g  .     ,  ■    j         r.     ■  •  en  ^T        9    inlet      f    inlet 

v/mdow  windov;   Regime  Samples      d£        ^ 


p      at 

g 

dndc 
loc .  loc. 


LB 
m 

HR 

psia 

LB 
m 

3 

FT 

psi 

FT 

°F 

°F 

14.42 

30.69 

0.14  8  6 

S 

0.606 

7 

0.527 

74.5 

63.5 

10.55 

29.93 

0.1429 

S 

0.596 

7 

0.4  87 

81 

61 

8.42 

29.86 

0.14  27 

S 

0.55  4 

10 

0.503 

81 

61 

6.22 

29.82 

0.1420 

S 

0.568 

6 

0.483 

83 

61 

4.075 

28.  28 

0.1375 

B 

0.504 

3 

0.450 

83 

61 

2.  690 

26.93 

0.1282 

B 

0.433 

3 

0.417 

83 

61 

1.490 

24.2 

0.1150 

B 

0.290 

3 

0.403 

83 

61 

1.470 

2  4.9 

0.1189 

B 

0.317 

3 

0.403 

83 

61 

1.018 

23.26 

0.1105 

B 

0.209 

3 

0.369 

83 

61 

0.  622 

22.46 

0.1071 

B 

0.146 

3 

0.376 

83 

61 

0.284 

21.32 

0.1015 

B 

0.056 

3 

0.413 

83 

61 

B-S    Transition    at   W      =    6.0    LB   /HR 

g         m 


Flow  Regime  Key 

S-A   Slug  annular 

S     Slug 

B--S   Bubbly-Slug  transition 

B     Bubbly 

10  4 


APPENDIX  D 
List  _' ?f  Equipment 

Pressure  Level  Measurement 

Stratham  Pressure  Transducer  Model  UC2 
Stratham  Model  UR5  Universal  Transducer  Readout 

Differential  Pressure  Measurement 

Stratham  11639  PM  60  tp  n>v.~       ~  •  ,- 

fm   bu  rc  Pressure  Differential  Transducer 

Ramapo  Strain  Gage  Adapter  Model  S6A  100B 

Hewlett  Packard  Model  A  i  •>*  t^  tt 

.uoaei  q]2A   DC  Vacuum  Tube  Voltmeter 

Air  and  Water  Flow  Measurement 

Fisher-Porter  Flowrator,  (various  models) 

Void  Fraction  Isolation  Valvos 

Rockwood  Figure  101  Bal]  Valves 

H    1/2"  -  main  valve;-  i  /a  "    u         -, 

5,1  -V4   ~  bypass  valve) 

Electric  Conductivity  Probe 

24  gauge  stainless  steel  spring  wire 
MIC-062-A  Conax  Fittings 

EPO-TEK  418H  -  „igh  Temperature,  Electrioally  Conductlve 

Epoxy,  (Kpoxy  Technology,  Inc.) 
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R277    Reese 

Void  fraction  in 
parallel  channel  boil- 
ing water  reactors  with 
application  to  two- 
phase  turbulent  mix- 
ing. 
24JUL70  nt^ptAY 
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R2?7         Reese 

Void  fraction  in 
parallel  channel  boil- 
ing water  reactors  with 
application  to  two- 
phase  turbulent  mix- 
ing. 
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Void  fraction  in  parallel  channel  boi'in 
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